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The basic principle of astronomical interferometry is to derive the angular distribution of radiation 
in the sky from the Fourier transform of the electric field on the ground. What is so special about the 
Fourier transform? Nothing, it turns out. I consider the possibility of performing other transforms 
on the electric field with digital technology. The Fractional Fourier Transform (FrFT) is useful for 
interpreting observations of sources that are close to the interferometer (in the atmosphere for radio 
interferometers). Essentially, applying the FrFT focuses the array somewhere nearer than infinity. 
Combined with the other Linear Canonical Transforms, any homogeneous linear optical system with 
thin elements can be instantiated. The time variation of the electric field can also be decomposed 
into other bases besides the Fourier modes, which is especially useful for dispersed transients or 
quick pulses. I discuss why the Fourier basis is so commonly used, and suggest it is partly because 
most astrophysical sources vary slowly in time. 


I. INTRODUCTION 

Aperture synthesis interferometers are now a common 
instrument in astronomy that map sources on small angu¬ 
lar scales. Interferometers exploit the time/phase delay 
of a radio signal arriving at different receivers to measure 
these details [I|. The typical understanding of how in¬ 
terferometers work is based on the Fourier transforms of 
the electric field - first in time, then in space 0. Inter¬ 
ferometers have been used for decades now, particularly 
in radio astronomy, where the long wavelengths lead to 
low angular resolution Q. 

Normally, diffraction limits the resolution of a single¬ 
dish telescope with diameter L observing radiation with 
wavelength A to A9 k, XjL. But a pair of receivers in an 
interferometer measure the Fourier mode of the angular 
distribution of radiation on the sky with characteristic 
angles 

Ad « A/£, (1) 

if £ is the separation of the receivers. Each measured 
mode - from each pair of receivers - corresponds to a 
point on the uv plane, which is basically the Fourier 
transform of the beam shape. A single dish telescope can 
be viewed as a filled-in array of receivers that measures 
all of the Fourier modes within a disk of the uv plane. 
While a pair of interferometers provides only sparse infor¬ 
mation, more detailed geometrical information can be re¬ 
constructed by using an array of N interferometers (with 
0{N^) pairs), using Earth’s rotation to rotate the uv 
plane, or using a spread of different frequencies. 

Although early interferometers actually channeled ra¬ 
diation from receivers to a central location (for example, 
through waveguides) to be correlated, with modern tech¬ 
nology, the field amplitude at each receiver can simply 
be recorded digitally. Then, the signals can be processed 
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numerically later to reconstruct the image. Digital pro¬ 
cessing of this kind is necessary for Very Long Baseline 
Interferometry (VLBI), where the receivers can be sep¬ 
arated by thousands of kilometers Q. It is the method 
of choice for new or upgraded radio telescopes, includ¬ 
ing the Karl G. Jansky Very Large Array (JVLA) [^, 
the Giant Metre-wave Radio Telescope (GMRT) Q, the 
Allen Telescope Array (ATA) Q , and the Low Erequency 
Array (LOFAR) @. Digital aperture synthesis is now en¬ 
abling several arrays that are designed to detect 21 cm 
emission from the high-redshift intergalactic medium Q , 
among which are the Donald C. Backer Precision A rray 
for Probing the Epoch of Reionization (PAPER) [lOj . 
the Murchison Widefield Array (MWA) 0 , the Aus¬ 
tralian SKA Pathfinder (ASKAP) [l^, and MEERKat 
[l^ . These are currently intended to be pathfinders for 
the proposed Square Kilometer Array (SKA), a global 
observatory intended to be built in the 2020s [TJ. 

The Fast Fourier Transform Telescope (FFTT; i) is 
a refinement of digital interferometers. The original pro¬ 
posal calls for a large N x N rectangular grid of simple 
dipole receivers {N = N'^ total) with spacing £ < A. 
The gridding of the antennas allows the electric field to 
be transformed with the Fast Fourier Transform (FFT) 
in only O(fVlogiV) operations S[I1. The FFTT can 
map radio emission on angular scales ranging from the 
whole sky to A/(£V) radians. It is thus ideal for measur¬ 
ing moderate to large brightness features over vast fields, 
the principle requirement for 21 cm reionization studies. 
The array of dipole antennas at the Maipii Radio Astron¬ 
omy Observatory was an analog version of the FFTT, al¬ 
though there were only 6 independent rows of antennas 
[T^ . The Waseda FFT Telescope is another early imple¬ 
mentation of the FFTT, with N = 16 [l^. The FFTT 
can be simplified into the fractal “omniscope” which is 
a nested series of grids while still providing resolutions 
over a wide range of scales [l^ . One of the low frequency 
pathfinders, the MIT Epoch of Reionization (MITEoR) 
telescope, is the first modern omniscope to be built [l8l| . 

But fundamentally, why do we use the Eourier trans¬ 
form in both time and space? Gosines and sines are 
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not the only possible basis for representing functions. In 
fact, the Fourier transform is only one of a whole fam¬ 
ily of transformations, the Linear Canonical Transforms 
(LCTs). These transforms do have physical meanings, 
which will be explored in this paper. They can perform 
any area- and orientation-preserving linear mapping on 
space-frequency phase space: rotations, squeezing, and 
shearing. Nor are we limited to LCTs: we could express 
the electric field as the linear sum of any set of basis 
functions, such as wavelets or pulses. In principle, even 
more general non-linear transforms could be performed 
digitally. 

Understanding the meaning of these other transforms 
provides insight into how interferometry and beamform¬ 
ing works. Although I have chosen to center my discus¬ 
sion on astronomical interferometers, the basic ideas ap¬ 
ply whenever beamforming is used to receive or broadcast 
a signal. Beamforming, of course, is the main principle 
behind phased array antennas, which are used in radar 
detection. Nor is beamforming limited to electromag¬ 
netic radiation, but it also is used frequently with sound 
waves, such as sonar. 

A note on notation: throughout this paper, I use S(x) 
to refer to the Dirac delta function, i refers to V—l, and 
an overline over a variable represents its complex conju¬ 
gate. 


II. HOW STANDARD INTERFEROMETRY 
WORKS: AN INTUITIVE EXPLANATION 


Interferometers exploit the phases of arriving electro¬ 
magnetic waves to reconstruct the direction the radiation 
is coming from. Consider polarized monochromatic elec¬ 
tromagnetic waves emanating from a single direction § 
on the sky.^ The electric field at a point r and time t is 


E = Eo exp 


27ri 

— (s • r -I- ct) 
A 


( 2 ) 


if the wavelength is A. I will consider the case of a ID 
array, but the basic ideas can easily be generalized to 
a 2D array. In addition, I will suppose that all of the 
receivers are located on the ground (2 = 0). The receivers 
in a ID array are located at different positions in the x- 
direction. The electric field on the ground is 


E(x) 


Eo exp 


2Tri , . „ , 

—— (xsint^ -I- ct) 

A 


(3) 


where 9 is the zenith angle of the source. 

If the source is at the zenith, the phase of the electric 
field is the same everywhere on the ground (Figure [TJ. 


However, as the zenith angle increases, it takes radiation 
longer to reach one side of the array than the other be¬ 
cause the path lengths are slightly different. The electric 
field on the ground at any instant therefore oscillates on 
the x-axis, with spatial wavelength = A/sin0 (Fig¬ 
ure [T|). The spatial Fourier transform of this electric field 
is particularly simple: 

J-[E](k)=EoX(5(^ac-^^ (4) 

In Fourier space, all of the power is concentrated at an 
inverse wavelength k with a simple dependence on 0.^ 
Thus, the Fourier transform of the held describes the 
direction the radiation is coming from (bottom panels 
of Figure [T]). Since the Fourier transform is a linear 
operator, additional sources sum linearly in the Fourier 
transform of the electric held. The interferometer can be 
“aimed” at a given direction Sq by adding a time delay 

At = §0 ■ r/c (5) 

to each receiver, adjusting the relative phases of the elec¬ 
tric held at each position. 

When considering radiation of many different wave¬ 
lengths, there appears to be an ambiguity in equation 2] 
How can we tell whether a short wavelength mode is a 
high frequency wave coming from near the horizon or a 
low frequency wave coming from near the zenith? This is 
easily resolved by examining the signal in time domain. 
The time frequency of the signal is the same everywhere 
on the ground plane; it is merely delayed at some points 
more than others. Interferometry works because the pro¬ 
jection of the wave onto a plane stretches it out in spatial 
wavelength without changing its time frequency. 

Each pair of receivers a and b is separated by a distance 
{Axab, Ayah)', the separation can also be written in terms 
of the uv coordinates: 

/ \ (AXab sin ^source Ay^b sin ^source 

\Uabi Vab) — I j ^ 

The pair measures a single Fourier mode of the electric 
field with wavelength (uo&A, UabA), which only provides 
information about the distribution of radiation on the 
sky on an angular scale of 

{9ui 9v')ab — (1/^ab; 1/^ab) ■ (?) 

Thus, traditional interferometers are insensi¬ 
tive to structures with angular scales greater 
than (min(l/uob), min(l/uob)) or less than 
(max(l/uab), max(l/uab)). An important goal of 



^ The actual sum of the sky radiation can be described as a linear 
sum of such waves with different frequencies and coming from 
different directions. 


^ I use K = 1/A instead of the wavenumber k = K/{2n) throughout 
the paper. This avoids the need to divide by 2it when considering 
the units, although the equations themselves then do use a factor 
of 2it. 
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FIG. 1: As the angular altitude of a source of plane waves increases, the projected spatial wavelength of the plane waves 
decreases due to projection effects. The components of the electric field are imaginary (red) and real (blue), with a magnitude 
shown iir black. The solid lines iirclude only the signal that lands on the array, while the dotted lines show how the signal 
continues beyond its bounds. 


interferometers is to “fill in” as much of the uv plane as 
possible, and this can be aided by allowing the Earth’s 
spin to rotate the sky and by using different frequencies. 

The FFTT is the simplest interferometer to under¬ 
stand without obscuring the basic principle. In a stan¬ 
dard FFTT, the receivers are arranged on a N x N rect¬ 
angular grid with separation Al < Amin- By Nyquist’s 
theorem, the electric field everywhere on the grid can be 
reconstructed with a spatial resolution Aniin/2. A Fast 
Fourier Transform can be directly be performed on the 
electric fields measured at each receiver to produce an 
electric field image, which is then squared to give the in¬ 
tensity image 0. Structures on angular scales ranging 
from Xmin/ (NX) to the whole sky can be mapped this 
way with the FFTT. 

The spatial Fourier Transform works because it mea¬ 
sures the projection effect on the time-domain Fourier 
modes. 


A. Complications in real interferometers 

Of course, it is often impractical to measure the electric 
field everywhere as a continuous function. The subtleties 
in understanding interferometry (particularly the use of 
“visibilities”) arise because of the limited sampling of the 
electric field. 

Actual interferometers are composed of a series of dis¬ 
crete elements, which could be anything from simple 
dipole antennas to single dish telescopes. These receivers 
each have a primary beam A[k) that describes their sen¬ 
sitivity to radiation that comes from different directions 
K on the sky. A large single dish’s primary beam is typi¬ 
cally confined to a narrow region in the sky with radius 
A/Hdish for a dish diameter D^ish- Dipole antennas have 
sensitivity over much of the sky, which is why they are 
proposed for FFTTs. The measured electric field (or in¬ 


tensity) is the sky intensity convolved with A{k) (or 
A{k)). I assume that the antennas all have the same re¬ 
sponse, are all pointed in the same direction, and there 
is only one point source in the sky, so that the antenna 
response can be approximated as a constant \/~A (or A 
in intensity). 

The other limitation, mainly faced by traditional in¬ 
terferometers, is that the uv plane may be sparsely sam¬ 
pled. Interferometers, whether FFTTs or traditional, are 
insensitive to the electric fields outside of their collective 
area. In addition, sparse interferometers only have a few 
measurements of the electric field from which to extrap¬ 
olate the entire uv plane. Thus, the measured electric 
fields are filtered by the interferometer’s sampling on the 
ground, S'(x); 

F;(i(x) = VAFn(x) X 5(x), (8) 

for a polarization 11. Then the Fourier transform of the 
electric field is the convolution of its Fourier modes with 
the dirty beam Pdirty = -Pis']: 

= F[E'n]{K) = VIP[£;n(x)] * . (9) 

The result is the dirty map. Some method is then neces¬ 
sary to deconvolve the dirty map with the dirty beam to 
get the sky distribution of radiation. 

In traditional interferometry, these difficulties are dealt 
with by calculating visibilities for each pair of receivers, 

V{ui2,Vl2) = {E 1 E 2 ), ( 10 ) 

averaged over time.^ The van Cittert-Zernike theorem 
relates the visibilities to the Fourier transform of the sky 


® Note that intensity (surface brightness) has units of the electric 
field squared. 
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intensity (e.g., M)- Specifically, the dirty map is calcu¬ 
lated as 


/dirty («) = AJ^[V] * -f^dirty- (11) 

Because of the sparsity of traditional interferometers, 
their dirty beams may be very complicated. Algorithms 
like CLEAN [2^ deconvolve the dirty map with the dirty 
beam. 


III. DIGITAL OPTICS THROUGH LINEAR 
CANONICAL TRANSFORMS 


A. A whole family of integral transforms 


The Fourier transform is actually just one of the family 
of linear canonical transforms (LCTs). The LCTs are the 
set of all linear integral transforms that conserve volume 
and orientation in position-frequency space. They can be 
represented by the ABCD matrix 


“=( cd ) 

with a determinant of 1 [Ulil. The LCT synthesizes 
a signal at a new position coordinate x' and a new fre¬ 
quency (inverse wavelength) coordinate k' out of the old 
position X and frequency k, as given by 



(13) 


There is no unique definition for an LCT, but the follow¬ 
ing is usually used [^ : 


\Mf{x)]{x') = 



J(x)dx 

(syo) 


(B = 0) 
(14) 

The Wigner distribution provides a geometric analogy 
for the effects of LCTs. It measures how much of the 
signal is at each position and frequency 


(15) 

If f{x) is a delta function, then its Wigner distribution 
is a vertical line (taking x to be horizontal and k as ver¬ 
tical). The function exp(27riKa;) is a horizontal line in 
Wigner space. Diagonal lines in Wigner space are linear 
chirps. 

A form of Parseval’s theorem applies to LCTs, 


nOO pOO 

/ |M/(x)p = / |/(x)p. (16) 

J (X> J oo 


As such, an LCT cannot increase the total “energy” in a 
signal. What they can do is carve up Wigner space into 
slices that are aligned with a signal’s energy. This con¬ 
centrates the energy into a small bin along one direction, 
making the signal easier to detect. 

The LCTs are simply the composition of a few ba¬ 
sic operations in Wigner (position-frequency) space (^ . 
The first basic operation is rotation through some angle 
a: 


Ja 


/ cos a sin a 
\ — sin a cos a 


(17) 


Rotation mixes the position domain with frequency do¬ 
main. When a = 7r/2, the LCT is simply the Fourier 
transform: 


T = 




(18) 


Then the new position coordinate is k. For other angles, 
the rotation is a Fractional Fourier Transform (FrFTs). 
FrFTs use chirps as their orthogonal basis functions. The 
second basic operation is shear by some quantity q: 


Fr — 

riq — 



(19) 


The shear of position/frequency space is also known as a 
Fresnel transform. The final basic operation is magnifi¬ 
cation by a factor fi, or 


= (o JJ ■ P") 

Nor are LCTs limited to real numbers. The Laplace 
transform is well known to be the Fourier transform of 
the imaginary frequencies, and is given by: 


C = 



( 21 ) 


Of course, other transforms can be constructed out of 
these basic operations. 


B. LCTs and matrix optics: The basic optical 
elements 

The LCTs are deeply connected with optics. Quadratic 
phase systems can be represented by an ABCD matrix 
that corresponds directly to the ABCD matrix of an 
LCT [ 2 ^. Suppose the electric field of a light wave at 
location 1 is described in terms of a position xi and spa¬ 
tial frequency ki. The light wave is processed by some 
optical system M and reaches location 2, where the elec¬ 
tric field is described in terms of position X 2 and spatial 
frequency K 2 . Then the optical system’s effects can be 
represented as 



( 22 ) 
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For monochromatic light, the spatial frequency directly 
corresponds to a ray’s angle with the optical axis. 

M is a composition of the ABCD matrices of several 
basic elements. If a light wave moves a distance d through 
free space, the matrix is a shear matrix 


M = 


1 d\ 
0 1 


(23) 


This is of course the Fresnel transform, named for Fresnel 
diffraction. A thin lens with focal length / shears in the 
frequency axis instead: 


M = 


:) 


(24) 


Any LCT can be constructed out of these elements [2^, 

With digital optics, there are none of the usual physi¬ 
cal limits on what material properties are realistic or not. 
One can easily emulate a material with index of refraction 
that is less than 1, negative, or even complex. Digital op¬ 
tics that mimic negative index of refraction materials can 
effectively act as superlenses, useful for near-held imaging 
of sources closer than one wavelength (see section IVT). 


C. Lenses as Fourier transformers 


As is well known, a lens performs the Fourier transform 
of the electric held. We can conhrm this with the ABCD 
formalism. Consider a lens with focal length /. It focuses 
light from inhnity into an image on a plane a distance 
/ after the lens. Take the input held to be located a 
distance / in front of the lens and the output held to be 
/ behind the lens. The ABCD matrix is 


l^focus 


1 A/ 

0 1 
0 

-1/(A/) 


1 0 
-1/(A/) 1 
A/ 

0 



(25) 

(26) 


This is just the Fourier transform LCT matrix done in 
units of i/A/. That is, we perform a “magnihcation” by 
a factor l/y/\J on the input plane to get a unitless co¬ 
ordinate system, perform a Fourier transform, and then 
perform a “magnihcation” by a factor i/A/ on the output 
to get it back into physical units. 


IV. FRACTIONAL FOURIER TRANSFORMS: 
HOW TO FOCUS AN INTERFEROMETER 



FIG. 2: The geometry of the fractional interferometer. 


the sky. The identity transform on the electric held is 
then most useful. If the source is at inhnity, it lies in a 
dehnite direction, but it has no “position” within the ar¬ 
ray. Then the hrst-order Fourier transform of the electric 
held is necessary to recover its direction. This thought 
experiment suggests that if a source is at intermediate 
distance, its position and direction are both moderately 
dehned, and one would use some intermediate fractional- 
order Fourier transform to hnd it. This turns out to be 
the case. 


A. Parallax and FrFTs 

If the source is merely very distant but not at inhnity, 
the source displays parallax between different parts of the 
interferometer (Figure [2|). Because of parallax, the radi¬ 
ation cannot be effectively described as coming from a 
single direction, even if the source is localized (Figure[2|). 
The spatial wavelength of the radiation is smeared out, 
and the Fourier transform normally implemented by in¬ 
terferometers performs poorly. 

This can be seen in detail by reconsidering the case 
of a ID array, but now with a point source located at a 
position X = —x^, and z = Zi, (Figure [2|). Its distance 
from the center of the array is V. The electric held at 
any given time is: 


Fractional Fourier Transforms lie in between the stan¬ 
dard hrst-order Fourier transform and a zeroth-order 
Fourier transform - the identity transform. One can 
think of the identity and Fourier transforms as two lim¬ 
iting extremes. If the source is so close that it lies within 
the array, it has a position but no single direction in 


E(a;,z,t) = 


X exp 


Eo 


Y^(x -k Xi,Y + {z- z*)2 
27rz , 


A 


■iy/{x + a:*)2 + {z - z*)2 -|- ct) 


(27) 


The electric held on the ground has a complicated de- 
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FIG. 3: The spatial frequency of waves from a nearby source 
changes with distance. The signal on the ground is then a 
chirp. This smears out the image in standard interferometers. 
The lines styles are the same as in Figure [T] 


pendence on position: 


E(x) = 


En 


exp 


/ 2'Kict\ 


^J{x + x^,Y +zl \ ^ J 


X exp 




(28) 


and thus cannot be inverted easily. The spatially oscil¬ 
lating part of E(a::) is an exponential of the form exp(i(f>). 
The phase (j) is 




2ttV 


1-h 


(X\ . „ X^ 

I — I sin 0 H- tt 


cos^ 6 


(29) 


when considering regions where x <^'D. 

Taking a spatial derivative of (j) amounts to measuring 
the instantaneous spatial frequency, 

27rK = ^ ^ ("sin0 -|- — cos^ 0 ) . (30) 

This frequency is what interferometers usually measure, 
and the Hrst term is just the projection effect described 
in Section m But if D is not inhnite, then the spatial 
frequency does depend on position and the electric field 
can be approximated as a linear chirp . The frequency 
slide is 


dn d^cj) 
^ dx dx'^ 


2-k 1 

Td 


cos^ 0 . 


(31) 


What we need is a transform that breaks the signal into 
chirps, measuring the power in a domain intermediate 



FIG. 4: A plot of the distribution of a chirp’s signal in Wigner 
space (thick black/blue line). An interferometer can focus on 
a source by rotating by an angle a so that the energy is all at 
one position x'. The coordinates are in natural units of A for 
X and 1/A for k. 


between position and frequency. The Fractional Fourier 
Transform [ttL[^. 

^a[fix)]{x) = Vl -icota X 

/ OO 

exp[7rI(cota — 2 esc a x'x + cota X^)]fix)dx- 

-00 

(32) 

is suited for this purpose. (I use x 3'iid x' h^re to 
emphasize that the position and direction variables are 
unitless.) The angle a can be interpreted as the order 
Q = 0 ( 2 / 71 ) of the FrFT. The identity operation has 
Q = 0 (a = 0), while the standard Fourier transform has 
Q = 1 {a = 7r/2). With my conventions, x' reduces to 
K for Q = 1. Although FrFTs are not (directly) used in 
astronomical interferometry to my knowledge, they were 
defined in 1980 with awareness sprea ding in fields such 
as optics (e.g., 1^1. signal processing [^, and quantum 
mechanics (e.g., [^) through the 1990s and 2000s. 

In Wigner space, the distribution of the energy in E(a;) 
is a diagonal line as shown in Figure |H The FrFT ro¬ 
tates position-frequency space so that the line of energy 
becomes vertical - a ^-function ISection llll All . Then one 
can just read off the chirp’s frequency derivative and use 
it to infer the distance of the source. But we must choose 
the “angle” a correctly in order to optimally concentrate 
the chirp’s power. 

The frequency k plotted in Figure |4] is l/(27r) x 
(d(j)/dx). Say the units for position are s and the units of 
K are then 1/s, in order to put everything into a unitless 
form (with x = x/s). The slope of the chirp in Wigner 
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space is 


tan 5 = 


d{Ks) 

d{x/s) 


d‘^4) 

27r 


cos^ 6 
XV 


(33) 


From this I infer that /3 = tan ^ s^/(AZ)) cos^ d], and 

the necessary rotation angle is 


a = 


^ ■S 2 /> 
TT + cos^ 9 . 


(34) 


In other words, we must perform a Fractional Fourier 
Transform of order Q Ri 1 + (2/7r)(scos0)^/(A2?). I call 
an interferometer that performs a FrFT a fractional in¬ 
terferometer. 

In astronomy, the sources are typically very distant 
and the units are fairly small, so that Q —>■ 1. The most 
“natural” units are the wavelength of the radiation A, 
in which case the deviation of Q from unity is a mere ^ 
A /V ^ I. Then x' is very nearly k and has practically the 
same units of inverse wavelength. In this limit, cot a —>■ 
—j5 + 0{fi^) and esc a —1 + The calculation of 

the FrFT is easy in this limit: 

J^a[E]ix') ^ expl-TTiPx''^] T[E{x) exp^-TTiPx^)], (35) 
or after restoring units, 


Ea[E]{x') « exp 


E 




COS 


X 

^2 


E{x) exp ( —ni^j^cos^O 


■ (36) 


A transfer matrix of this type can be converted into a 
FrFT matrix by a scaling with the correct units: 

/ I /1//7 ® ^ M sin a A / ^Tfi 0 

“ \ 0 y/7/ \ —(1 /m) sina cosa y V 6 

“ (38) 

In this case, the effective angle of the FrFT performed 
by the optical system is cos a = 1 — d/ f if we put 
position in units of ^/Xf (compare section ITlI Cl) . 

When a source is paraxial with the lens (cos 9=1) and 
is located a distance V in front of it, its image is found 
a distance 


d = f/{l-f/V)^f{l + f/V) (39) 

behind the lens. The imaging process is then equivalent 
to a FrFT with order given by cos a ~ —fjV. We can 
compare this with the FrFT order given by equation 1341 
In that case, if the unit s is chosen to be v^A/, then 
cos a = — sin/3 Ri —ffV, an identical result. 

An interferometer that performs the FrFT on an elec¬ 
tric field focuses an image. Conventional interferometers 
are focused at infinity, as suited for most astronomical 
sources. But there is no reason why they cannot focus 
on things that are closer. 


C. Fractional dirty beams, fractional convolutions, 
and related horrors 


These approximations apply under the loose condition 
that P'^xx' 1- Since the largest x in the problem is 
generally the array size A, and the largest x' is related to 
the minimum separation of the antennas {x' < I/^), this 
is true if A f(AA>)^/s^. If one works in a “natural” unit 
system suggested above, with s = A, then the condition 
is almost trivially true, A i{D/X)^. 

Applying a FrFT instead of a FFT therefore basically 
amounts to a chirp multiplication before and after the 
Fourier transform. 


B. The Fractional Fourier Transform as focus 

Just as a lens performs an ordinary Fourier transform 
when it focuses light from infinity (section fill Cl) , it per¬ 
forms a FrFT when it focuses light from distance V. 
Specifically, if we again consider a lens with focal length 
/, then the electric field in an output plane a distance d 
behind the lens is a Fractional Fourier transform of the 
field in an input plane a distance d in front of the lens 
[S^- This can be seen with a matrix optics analysis. The 
transfer matrix for light propagating distance d, passing 
through a lens with focal length /, and propagating d 
again is [U, [2^ 

_ /1-d// dA(2-d//)A 

'o™^“Ul/(A/) l-d/f )■ 


As with regular interferometers (Section fn]), measure¬ 
ments become much more complicated with actual inter¬ 
ferometers. First, the measured signal at a given antenna 
depends on its primary beam’s sensitivity in the direc¬ 
tion of the source. Frequently, the primary beams can 
be approximated as the same for each antenna. But for 
a nearby source, the parallax effect implies that the sen¬ 
sitivity of some antennas may be greater than others, if 
only because the antennas are not all pointed directly at 
the source. This could complicate the analysis if the par¬ 
allax angle is comparable to or larger than the primary 
beams. For distant objects with small parallax angles de¬ 
tected by dipole antennas, this may not be as much of an 
issue. As in section ITl A[ I will treat the antenna response 
as a constant -x/A, although this needs to be treated for 
more realistic situations. 

But a fundamental difficulty is the “filtering” 5'(x) per¬ 
formed by the interferometer as it samples the electric 
signal. The measured electric field is multiplied by S'(x), 
as in equation ISl and the Fractional Fourier Transformed 
field is 

Adirty = E4E'{x)] = E4VAE{x) X ,5(x)] (40) 

The usual convolution theorem of the Fourier transform 
does not apply (after all, the identity operation is Eq). 
Instead, one must replace the convolution operation with 














a fractional convolution operation that merely reduces to 
normal convolution when a = tt 12. Then we can say that 

i^dirty = = Pa[P(x)] pd"*^ 

(41) 

where p^irty _ r jg fractional dirty beam of the 
interferometer |23l| . Even if S' is a simple rectangular 
pulse, the fractional dirty beam does not have a sim- 
|?1^ analytic expression and it has imaginary components 

One example of a definition for fractional convolution 
is given by 

(/*o5)(x) = 

(42) 

which I have altered to fit this paper’s conventions re¬ 
garding factors of 27r (see also [^). Then, one can write 

Pdirtye”^'™*“ = (VAe-'^'™*“P«[P(x)]) 

* (e’^^x'cotapdirty)^ (43) 

which is a standard deconvolution problem. 


But if the source is too close to the array, then the 
approximation in equation [221 is insufficient. The next 
term in (j) is {2'k'D/X) x [—(l/2)(x/P)^ sindcos^ 0]. The 
quadratic approximation is sufficient only if this next 
term adds or subtracts less than 1 radian to the phase: 


sm0cos^0 


> 1 . 


(46) 


We have 


V < L\ cos( 


7 rL| sindj 


A 


= 32 km I cosdj i/j sindj 



The above condition is listed for various arrays in Table U 
too, under ^quadratic- Even higher mth order terms of 
the approximation are proportional to a;™/(AP"*“^), and 
are valid for V < ^ ^j^g Xaylor series 

breaks down entirely when 7) < L. 


E. Subarrays or fractional interferometry? 


D. The limits of fractional interferometers 


Unless one works at very low frequencies (e.g., study¬ 
ing atmospheric phenomena at kHz frequencies), we are 
typically in the regime where P 3> A. This condition 
motivates my use of the Taylor series for P'/P in equa¬ 
tion [221 The reader may notice that the order Q of the 
FrFT is very close to 1 in this case 1 equation (341) . Al¬ 
though one might worry whether the interferometer can 
measure the phase of the electric field well enough for 
the FrFT to be needed, it turns out they can for sources 
within the atmosphere and into cislunar space. 

An interferometer can measure the parallax effect of a 
nearby source if the parallax angle between one end of the 
array and the other is bigger than the limiting resolution 
from Rayleigh’s criterion. Another way to think about 
it is that the frequency derivative of the chirp can be 
measured if the phase slides more than 1 radian across 
the array from a constant tone. This gives us: 


p Ri tanp 


LcosO ^ A 
P L 


(44) 


for an array of width L. The relevant condition is 


T* < 


cos 9 
A 


= 330 km cos 9 


km.J \ 100 MHz 


(45) 

using values that could be typical of a 21 cm survey. 
When using the longest baselines of LOFAR or SKA (^ 
1000 km), or the highest frequencies of VLA (~ 50 GHz), 
the Moon itself is just out of focus with a blur of ~ 0.6 
arcsec (about 1 km). These limits for various arrays are 
listed in Table 11 under Punear- 


Fractional interferometry corresponds to a second- 
order approximation of the phase with position, or a first- 
order approximation in spatial frequency. This seems ob¬ 
viously advantageous over normal interferometry, which 
is a first-order approximation of the phase (zeroth-order 
in frequency) - we fit the signal with a line of arbitrary 
slope instead of a zero-slope line in Wigner space (Fig¬ 
ure a. But another way to describe a signal is to break 
the array into many smaller arrays and fit the signal in 
Wigner space with a series of constant values, one for 
each sub-array. Likewise, one could try to fit a signal 
with a nonzero second derivative of k by breaking the 
signal into a series of lines, to get around the limits im¬ 
posed by equation im 

Essentially, the problem is to minimize the angular 
beam width W if the source may be out of focus. Con¬ 
sider a ID subarray centered at x = xq with a total width 
of 1. There are two terms in W. The first is the reso¬ 
lution limit imposed by diffraction, Wdifi = A/Z, which 
favors large £. The second is the apparent angular size 
of the source’s image, Wimage, which gets larger as the 
source goes out of focus. Wimage is basically the projected 
“width” (or overhang) of the signal in Wigner space after 
the interferometer performs a FFT or FrFT; it is calcu¬ 
lated as the overhang in x' (which is almost, but not 
technically, k) in units of 1/A. With a standard inter¬ 
ferometer, the best-fit linear approximation to the phase 
leaves a residual signal width of 


TTylinear 

image 


l cos^ 9 
2V 


( 3 ^ sin 6» -k 1 ) 


I cos^ 9 
2V 


(48) 


for a distant source. (Note that Xq ranges from 
—1/2 to 1/2.) Minimizing this term favors small 1. 













9 


Array 

L 

A 

V 

^linear 

^quadratic 

References 

Kilometric radio 

Astronomical Low Frequency Array (ALFA) 100 km 

10 km 

30 kHz 

1000 km 

560 km 

[33J 



10 m 

30 MHz 

10® km 

1.8 X 10'‘ km 


Hobart array 

1 km 

60 m 

5 MHz 

17 km 

7.2 km 

[341 

Meter radio 

LOFAR Superterp 

240 m 

10 m 

30 MHz 

5.8 km 

2.1 km 

[8] 



1.5 m 

200 MHz 

9.6 km 

1.9 km 


LOFAR Core 

3.5 km 

10 m 

30 MHz 

1200 km 

120 km 




1.5 m 

200 MHz 

8200 km 

300 km 


LOFAR Remote 

121 km 

10 m 

30 MHz 

1.5 X 10® km 

2.4 X 10"* km 

[8] 



1.5 m 

200 MHz 

0.065 AU 

6.1 X lO'^ km 


LOFAR All 

1158 km 

10 m 

30 MHz 

0.90 AU 

7.0 X 10® km 

[8] 



1.5 m 

200 MHz 

6.0 AU 

0.012 AU 


Centimeter radio 

GBT 

110 m 

300 cm 

100 MHz 

4.1 km 

1.2 km 

mi 



0.26 cm 

115 GHz 

4700 km 

40 km 


VLA (D configuration) 

1 km 

30 cm 

1 GHz 

3300 km 

100 km 

m 



0.6 cm 

50 GHz 

1.7 X 10® km 

720 km 


VLA (A configuration) 

36 km 

30 cm 

1 GHz 

0.029 AU 

2.2 X lO'^ km 

mi 



0.6 cm 

50 GHz 

1.4 AU 

1.6 X 10® km 


VLBl 

10'‘ km 

20 cm 

1.5 GHz 

0.016 pc 

0.84 AU 




1 cm 

30 GHz 

0.32 pc 

3.7 AU 


RadioAstron -|- VLBI 

3.5 X 10® 

km 92.5 cm 

324 MHz 

4.3 pc 

81 AU 

[37] 



1.35 cm 

22.2 GHz 

290 pc 

670 AU 


Interplanetary radio interferometer 

1 AU 

20 cm 

1.5 GHz 

3.6 Mpc 

7.4 pc 


Submillimeter 

Event Horizon Telescope 

10'‘ km 

1 mm 

300 GHz 

3.2 pc 

12 AU 

[381 

Near infrared and optical 

Thirty meter telescope 

30 m 

0.6 /rm 

500 THz 

0.010 AU 

380 km 


CHARA 

331 m 

0.6 /rm 

500 THz 

1.2 AU 

1.4 X 10^ km 

[39] 

NIR VLBl 

10"^ km 

1 /rm 

300 THz 

3.2 kpc 

370 AU 


X-rays 

Maxim Pathfinder 

1.4 m 

200 pm 

150 PHz 

6.6 AU 

2100 km 

[40] 

Maxim 

300 m 

200 pm 

150 PHz 

1.5 pc 

0.044 AU 

m 


TABLE I: Closest distances for which normal and fractional interferometry apply for various arrays. These correspond to the 
domains where linear (< Ounear) and quadratic (< "Pquadratic) approximations to the phase with distance (equation 1291) are 
valid. 


Then, supposing the total beam width is about W ~ 
min[Wdiff, the optimal beam width is, 


TT/linear 

^^best 



(49) 


leaves a residual image width of 




3 

8 



sindcos^ 9. 


This allows an optimal width of 


for an array size of 

Ihesr ~ sec 9\/2XD. (50) 


w: 


quadratic 


best 



1/3 


Using a larger subarray (or telescope) than actu¬ 

ally makes the image blurrier because the parallax of the 
source across the array is detectable. 

The best-fit quadratic approximation to the phase 


for an array size of 


^quadratic 

best 


/o \l/3 

( csc9sec^ 9j 


(51) 


(52) 


( 53 ) 
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Again, a larger subarray than just makes things 

worse. 

Using a FrFT instead of a Fourier transform improves 
the image angular width as 


w: 


quadratic 


best 


TT/linear 

best 



1/6 


(54) 


for subarrays with the optimal size, where 


^quadratic 
best _ 

/linear 

^best 



1/6 


(55) 


Although these ratios evolve slowly with distance, note 
that the relevant factor of V/A is very big. As a result, 
the fractional interferometer can improve the best possi¬ 
ble beam width by orders of magnitude. 

I demonstrate this by calculating the best possible 
image resolutions with standard and fractional inter¬ 
ferometry for various sources in Table m For sources 
within the Solar System, standard interferometers have 
blurred, sub-optimal images for VLBI baselines. In addi¬ 
tion, if a RadioAstron-like baseline is ever added to the 
EHT, standard interferometry cannot resolve scales be¬ 
low ~ 1% of the Galactic Center’s supermassive black 
hole Schwarzschild radius. The fractional Fourier inter¬ 
ferometer performs much better. One would need base¬ 
lines that are spread across the Solar System before one 
would need to go to third order, if the sources are inter¬ 
stellar. Fractional interferometers are sufficient for our 
present astronomical needs. 


F. An example calculation of a fractional 
interferometer 

I shall now demonstrate how much using a FrFT in¬ 
stead of a simple Fourier transform matters for the beam 
size. Consider a ID FFTT that is L = 1 km wide de¬ 
signed to observe GHz radio waves. It has 2^® = 32768 
elements, spaced 3.05 cm apart. It therefore has enough 
information to sample I GHz radiation (A = 29.98 cm). 
I consider a case without any noise for simplicity. 

Now let’s say that one wishes to focus on a I GHz 
emitter that isV = 500 km away at a zenith angle of 45°. 
The transmitter may be a low earth orbit satellite, like 
the International Space Station (ISS), but it is a point 
source. According to equation [Ml the optimal FrFT has 
order very near I, with AQ = Q — 1 = 1.91 x 10“^. How 
much of a difference does this make? 

We can measure the effective image size by considering 
the smallest range of x' around the point of maximum 
intensity that contains 50% of the power. This quantity, 
dubbed IU 1 / 2 , is plotted in Figure [5] There is a flat 
bottom around the optimal AQ, but the image width 
begins increasing quickly as the deviation from this best 
AQ passes (0.5-1) x 10“^. A simple Fourier transform 
double the image size, and cuts the peak intensity in half. 


Note how the intensity peak is smeared out by the normal 
Fourier transform compared with the FrFT in the right 
panel. 

The fractional interferometer clearly performs better 
than the normal FFTT analysis, doubling the signal even 
for this tiny of a AQ. 


G. Targets for a fractional interferometer 


A fractional interferometer can increase an interferom¬ 
eter’s sensitivity and ability to maps sources that are rel¬ 
atively nearby. The inability of conventional interferom¬ 
eters to “focus” on nearby sources has been understood 
for a long time. Various fixes for focusing on these objects 
have been proposed, including adding a quadratic time 
delay for different elements (j^; see als o [^ i. and us¬ 
ing a nonlinear mapping on the uv plane [43| . The FrFT 
provides a natural framework for solving the problem. 

As seen in Table HI the entire Solar System is out of 
focus in radio VLBI arrays. Yet VLBI is an important 
tool for tracking space probes within the Solar System. 
Aside from simply keeping tabs on where the probes are, 
precise measurements of space probes are useful for pin¬ 
ning down the locations of the planets [4^, performing 
tests of general relativity and modified gravity (4^ . and 
monitoring the winds of planets and moons with atmo¬ 
spheric probes [H, . Among the oWects studied with 

VLBI tracking are Venus [i^. Mars [ij, Saturn Q, and 
Titan [13. 

The Moon is also a frequent target of VLBI obser¬ 
vations. The later Apollo missions planted radio bea¬ 
cons on the Moon’s surface as part of the Apollo Lu¬ 
nar Surface Experiments Packages (ALSEPs). Early 
selenometry benefited from VLBI observations of these 
beacons, which pinned down their relative positions and 
constrained the Moon’s motion VLBI tracking of 

lunar satellites like SELENE allowed for the determina¬ 
tion of the Moon’s gravitational field to high order • 

Artificial satellites in geosynchronous and low earth 
orbit are out of focus in normal VLBI observations, as 
well as in high-frequency observations by the Very Large 
Array (Table U) . These artificial satellites are monitored 
on occasion with long baseline radio interferometry [^ . 
There has also been recent interest in using optical inter¬ 
ferometry to directly observe communications satellites 
in geosynchronous orbit, to diagnose their status [^ . 
Yet, according to Table HI optical observations of these 
satellites are out of focus, and FrFT techniques must be 
applied to get a sharper image. 

Even closer to home, there are many radio phenomena 
in the atmosphere that could benefit from fractional in- 


The ALSEPs were turned off in the 1970s, but there’s no reason 
why future crewed or uncrewed missions to the Moon could not 
place additional radio transmitters. 
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Source 

V 

A 

Standard interferometer 

Fractional interferometer 




ilinear 

^best 

TTT-linear 

^^best 

A c*linear 
^•^best 

Tquadratic 

best 

Ti rquadratic 
^best 

^ ^quadratic 

Thunderstorm 

10 km 

10 m 

0.4 km 

1° 

0.2 km 

1 km 

O 

o 

70 m 

ISS 

500 km 

30 cm 

0.5 km 

2 arcmin 

0.3 km 

6 km 

10 arcsec 

30 m 

Geosynchronous satellite 3.6 x 10^ km 

1 ^tm 

8 m 

20 mas 

4 m 

2 km 

0.1 mas 

2 cm 

Moon 

3.84 X 10® km 

1 m 

30 km 

7 arcsec 

10 km 

700 km 

0.3 arcsec 

0.5 km 



1 nm 

30 m 

7 mas 

10 m 

7 km 

30 /ras 

5 cm 

Neptune 

30 AU 

1 m 

3000 km 

70 mas 

2000 km 

4 X 10® km 

0.5 mas 

10 km 



1 nm 

3 km 

70 /ras 

2 km 

4000 km 

60 nas 

1 m 

Solar neighborhood 

10 pc 

1 m 

8 X 10® km 

0.3 mas 

4 X 10® km 

4 AU 

0.3 /ras 

500 km 



1 fim 

800 km 

0.3 /ras 

400 km 

0.04 AU 

30 pas 

50 m 

Galactic Center 

8 kpc 

1 m 

0.1 AU 

9 /ras 

0.07 AU 

400 AU 

4 nas 

5000 km 



1 mm 

7 X 10® km 

0.3 /ras 

4 X 10® km 

40 AU 

40 pas 

50 km 



1 fim 

2 X 10"^ km 

9 nas 

1 X 10"* km 

4 AU 

0.4 pas 

0.5 km 



20 nm 

3000 km 

1 nas 

2000 km 

1 AU 

30 fas 

30 m 

M87 

20 Mpc 

1 m 

7 AU 

0.2 /ras 

4 AU 

0.3 pc 

20 pas 

6 X 10'‘ km 



1 mm 

0.2 AU 

6 nas 

0.1 AU 

0.03 pc 

0.2 pas 

600 km 



1 nm 

1 X 10® km 

0.2 nas 

5 X 10® km 

700 AU 

2 fas 

6 km 



20 nm 

2 X 10® km 

30 pas 

8 X 10"^ km 

200 AU 

0.2 fas 

0.5 km 

Distant galaxies 

2 Gpc 

1 m 

70 AU 

20 nas 

40 AU 

7 pc 

1 pas 

3 X 10® km 



1 mm 

2 AU 

0.6 nas 

1 AU 

0.7 pc 

10 fas 

3000 km 



1 /rm 

0.07 AU 

20 pas 

0.04 AU 

0.07 pc 

0.1 fas 

30 km 


TABLE II: Best possible angular resolutions with a standard interferometer (first order in phase) and a fractional interferometer 
(second order in phase). The largest resolvable features on the source is denoted by AS. All trigonometric terms on 0 are 
assumed to be ~ 1. 




AQ X 1 0^ A0 (arcmin) 

FIG. 5: How the characteristics of the image of the transmitter vary with FrFT order AQ (left), including the half-power beam 
width and the maximum intensity. On right, the intensity distributions of the image; plotted images are for AQ = 1.9 x 10~^ 
(solid), AQ = 0 (dashed), and AQ = —1.9 x 10~^ (dotted). The signal is normalized so that the total detected power is 1. 
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terferometry. Meteors vaporize into plasma in the atmo¬ 
sphere, which reflects artificial radio signal s . Meteors 
can emit radio emission on their own, too |54l | . Lightning 
frequently is a source of radio emission from Extremely 
Low Frequencies (< kHz) up through hundreds of MHz 
[ 5 ^ [ 5 ^. Sprites, which accompany some thunderstorms, 
also are sources of Very Low Frequency emission [^ . 
These kinds of electrical phenomena excite sferics and 
whistlers in the ionosphere at frequencies below a MHz 
[53 • Finally, Ultra High Energy Cosmic Rays (UHE- 
CRs) produce particle showers when they interact with 
the atmosphere. The detection of radio pulses accompa¬ 
nying these showers is an important technique for mea¬ 
suring the UHECR flux at the highest energies [H, . 

The high frequency radio transients in the atmosphere 
are of interest to the numerous radio astronomy arrays, 
and interferometry is occasionally used to study the low 
frequency transients associated with lightning . Frac¬ 
tional interferometers of the right frequency range may 
be more sensitive to all of these phenomena and could 
potentially map the shapes of the radio-emitting regions, 
perhaps providing new diagnostics. 

Observations of the new class(es) of Fast Radio Bursts 
(FRBs) are a good example of what a fractional interfer¬ 
ometer might accomplish. These millisecond-long pulses 
observed at GHz frequencies are dispersed in a way con¬ 
sistent with being located outside of the Milky Way [blj . 
But the situation is confused by the recent discovery of 
perytons, atmospheric radio transients of similar disper¬ 
sion, frequency, and duration . If a FFTT could catch 
some of these pulses, it could discern whether the FRBs 
are really astronomical or if they are atmospheric easily. 
An atmospheric FRB is out of focus [i^, and its image 
would be sharpened by a FrFT. Extraterrestrial FRBs, 
on the other hand, are essentially at infinity. 

Perhaps the most prosaic, but useful, application of the 
fractional interferometry technique would be to remove 
radio frequency interference (RFI) from human activities 
RFI is generated by electronic devices, both 
on the ground and in orbit. RFI is an ever-increasing 
threat to radio astronomy, and must be removed from 
radio data [i^ . In many cases listed in Table HI namely 
those with long baselines or at high frequency, the RFI 
sources are close enough to be out of focus. Focusing 
an interferometer could produce a sharper image of the 
RFI source, allowing it be removed more cleanly from the 
data without losing the wanted signal. 


V. THE USES OF EVANESCENT WAVES 

In the last section, I assumed that the radiation is high 
frequency, with both the array and the distance to the 
source much bigger than the wavelength. But at very low 
frequencies, it is quite possible that the array and the 
source distance are smaller than one wavelength. This 
could happen when observing at kHz frequencies, where 
the wavelength is 30 km (j^/IO kHz)~^, a common regime 


for sferic studies. Then we need to consider a near-held 
limit of microscopy. The recent advances in superlenses 
are applicable to this regime. 


A. Digital superlenses: perfect radio microscopes 

As I noted in Section llH Bl digital optics are not bound 
by normal physical constraints. Digital telescopes can 
simulate elements with negative index of refraction as 
easily as elements with positive index of refraction. Neg¬ 
ative index of refraction materials (NIMs) are especially 
interesting from an optical point of view, as they can be 
used to construct super lenses. 

Superlenses are lenses that can image objects that are 
smaller than the wavelength of the light they scatter as 
long as they are in the near-held limit (less than one 
wavelength away) Ira. Even a simple slab of NIM can 
act as a superlens [63 • When light enters a NIM, it bends 
away from the surface normal, according to Snell’s Law. 
This allows them to focus light both inside the material 
and again after it leaves the NIM slab [h^. 

The overall effect of a NIM is to effectively undo the 
propagation of light [h^. The reverse propagation ap¬ 
plies not only to the typical oscillating waves of light, 
but also exponentially decaying evanescent waves. The 
evanescent waves contain information about their source 
on scales below the wavelength. They vanish in the far- 
held {d> X), but if they enter a NIM before that, they are 
amplihed back up to their original strength or beyond. 
Thus, it is possible to reconstruct the electric held at the 
source in precise detail with a NIM. This superlens effect 
allows an optical system to capture information below 
the diffraction limit in the near-held limit [^, . 

Physical microwave superlenses have been built out of 
metamaterials, but there are fundamental limits to their 
performance that do not apply to digital optics. The 
most important is causality. A wavefront cannot bend 
away from the surface normal instantaneously; the far 
side of the wavefront needs time to catch up to the near 
side. Light in a NIM has a superluminal phase veloc¬ 
ity, but the group velocity is necessarily less than c, and 
there is a delay time before the wavefronts are launched 
0. The constraint of causality necessarily introduces 
dispersion, however, so a physical superlens is limited to 
a narrow wavelength range [68l| . 

But a digital system is not limited by causality. When 
a digital optical system is synthesized and applied to 
stored data, the entire history of the electric field is al¬ 
ready known, so an acausal system is possible. Causality 
is respected in a fundamental sense, because it already 
takes a light-crossing time for the information to be re¬ 
ceived by a central processor. As such, a digital superlens 
may be designed to have no dispersion and can work at 
any frequency. However, the resolution of a digital super¬ 
lens is still limited by the sampling of the electric field 

Digital superlenses are practical in the near-held 
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regime, when the source is less than one wavelength away. 
While this does not apply for astronomical sources, it 
does apply to atmospheric radio phenomena at very low 
frequencies. 

Astronomical sources are, of course, so far away that 
evanescent waves have long since decayed away. Trying 
to amplify these vanished waves instead merely amplifies 
random noise. The one proposal to use NIMs in astro¬ 
nomical telescopes involves artihcially creating evanes¬ 
cent waves in a NIM that simulate the evanescent waves 
that exist in the near-held of a source, which usually re¬ 
quires some knowledge of what the source looks like [t^ . 
An interferometer with a digital superlens is really a mi¬ 
croscope rather than a telescope. 


B. Laplacian interferometers: the (practically) 
impossible dream 


If digital superlenses use evanescent waves to recon¬ 
struct sub-diffraction images of near-held objects, it is 
natural to wonder if the evanescent waves can be detected 
directly. The evanescent waves are exponentially damped 
with distance from the source. The array measures a slice 
of these helds, which generically has an exponential de¬ 
pendence on position due to projection (Figure [5]). This 
is conceptually similar to normal interferometers, except 
that the exponentials of imaginary frequency are replaced 
by exponentials of real frequency. The necessary trans¬ 
form is the Fourier transform at imaginary frequency - 
the inverse Laplace transform. Can we build a Laplacian 
interferometer around this principle, instead of the usual 
Fourier interferometer? 

There are two types of Laplace transform that we can 
consider. The LCT formulation of the Laplace transform, 
given by 


£2 = 



(56) 


actually is the bilateral Laplace transform (aside from a 
phase factor): 

/ OO 

(57) 

-OO 

We also have the LCT version of the inverse bilateral 
Laplace transform: 

/ OO 

(58) 

-OO 

A single evanescent wave has a single projected scale kq, 
and has an electric field represented as a spike S{k— kq). 
The spike is Laplace transformed into an exponential in 
s-space, exp(—27rsKo). Ideally, we would just invert the 
measured electric field back into a delta function in n- 
space. Yet the inversion fails because the inverse bilateral 
Laplace integral diverges. Bilateral Laplace transforms 
generally have limited domains of convergence in C. 


The one-sided Laplace transform, 

nOO 

F{s) = £[f{K)] = / e-^y{K)dk (59) 
Jo 

has better convergence properties. Once again, the 
Laplace transform of 5{k — kq) is an exponential in s- 
space, exp(—s/to): as long as kq > 0. An exponen¬ 
tial can be inverse Laplace transformed back into a S- 
function. Conceivably, all we need to do to build an 
evanescent wave interferometer is perform a one-sided in¬ 
verse Laplace transform on the measured electric field. 

This is the infamous problem of separating exponen¬ 
tials that crops up from time to time in various fields 
[ 7 ^ . The underlying obstacle is that that the inversion 
requires integration in the complex s domain: 

r-i+iC 

/(/c) = £~^[/(k)] = lim / e^'^F{s)ds. (60) 

The contour of integration is the Bromwich contour, for 
which 7 has a greater real value than all singularities 
of F{s). While this is fine if we know the behavior of 
F(s) at complex s (as we do for the analytic exponential 
function), an array can measure the electric field only at 
real-valued positions. 

Because of the importance of the problem, several nu¬ 
merical algorithms have been developed that attempt to 
fit exponential(s) to data at real-valued positions [t^. 
For example, one might fit the data with an analytic func¬ 
tion (such as Laguerre or Chebyshev polynomials) and 
then extrapolate it into the complex domain [^ . The 
exponentials can also be fit by transforming the variables 
and performing a certain deconvolution [75| . 

All techniques suffer from extreme instabilities in the 
presence of noise, because the separation of exp onentials 
is “ill-posed” with several solutions [zlHiTg. In func¬ 
tional analysis terms, the exponential functions are far 
from orthogonal Since exponentials vary rapidly, 

even a slight fluctuation requires huge changes in the 
spectrum in order for the exponentials to cancel out in 
most places. Even the number of exponentials that pro¬ 
duce the signal can be impossible to extract [t^. 

The inability to infer which exponentials contribute to 
the electric field is fundamentally related to the diffrac¬ 
tion limit in optics , and why we cannot detect evanes¬ 
cent waves from astronomical sources. An image of an 
object is formed by an optical system with some entrance 
pupil. The image is a Fourier transform of the electric 
field near the pupil; therefore, the finite size of the pupil 
filters out high spatial frequency information. But the 
Fourier transform necessarily produces an analytic func¬ 
tion, so all of the information can hypothetically be re¬ 
constructed from the small sliver that passes the pupil 
[^ . The apparently infinite resolution compared to the 
limited resolution of actual optics presents an apparent 
paradox [zilzl- 

The paradox is solved by noting that although the high 
frequency information is still mathematically present, it 
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FIG. 6: Projection effects apply to a static electric field that has exponential scaling towards some direction in the sky. 
Separating these kinds of exponentials is impossible in practice. The lines styles are the same as in Figure [T] 


is suppressed to nearly zero amplitude. In more formal 
terms, the eigenvalues of the optical system’s eigenfunc¬ 
tions drop from almost exactly one to almost exactly zero 
as one passes the diffraction bound [zl, [t^ (see for 
a review). Attempting to find the strengths of the high 
frequency information results in a division by (nearly) 
zero, and even the slightest noise results in wild devia¬ 
tions. The problem is a generic one for integral equations, 
including the Laplace transform, which smooths the ex¬ 
ponential spectrum on small scales lzi,|8li. 

Physically, the suppression of the high eigenfunctions 
comes from the decay of evanescent waves decay after 
propagating any significant length. In order to measure 
those waves, they have to be multiplied back up by an 
exponential factor to their original strength - or, equiv¬ 
alently, their amplitudes must be divided by nearly zero. 
The inevitable noise within the system is also amplified 
with the evanescent waves, scrambling any information 
they contain. Even if we had a perfect detector, quantum 
fluctuations inject unavoidable noise, limiting the resolu¬ 
tion in the far-field [s^. Likewise, the throughput of the 
optical system is very small for evanescent waves, and 
any detectable signal requires a vast number of photons 
with a vast amount of energy [s^. Or, reversing time, an 
antenna can broadcast radiation in an arbitrarily narrow 
angle only at the cost of vast input power [sj . 

For astronomical objects, with distances® of 10®A, 
the evanescent waves are suppressed by a factor » . 

Hence, we cannot use super lenses to image astronomi¬ 
cal objects with arbitrary resolution. Only in the near- 
held, where the exponential functions are well-behaved 
and the evanescent waves have not yet decayed, can the 
Laplace transform be inverted. The virtual impossibil¬ 
ity of a Laplacian interferometer makes the connection 


® The Moon at 1 MHz. 


between the two problems explicit. 


VI. BEYOND THE FOURIER BASIS IN TIME 
DOMAIN 


A. Generic waveforms in the time domain 


Just as non-Fourier transforms can be used for the spa¬ 
tial structure of the electric held, the temporal structure 
of the held can be decomposed using other transforms. Of 
course, there is nothing inherently wrong with using the 
Fourier transform, as any function can be stably repre¬ 
sented as a sum of Fourier modes. The Fourier transform 
is not as natural, though, for representing signals whose 
frequency changes rapidly in time (as noted in [s^. 1^1. 

Suppose the radiation from a source is isotropic. We 
denote the electric held at a distance V (see Figure [5]) 
and a time t as E{T>'^t). Because the radiation is 
isotropic. 


E{V',t) 


m 

v 


f{t-V'/c) 

v 


(61) 


where f{t') can be any function with hnite energy. Now 
we wish to understand how the / waveform varies across 
the array. Let to = t — T)jc. Then, 


V 

t -= to 

c 


(P-P').,. a. 

-to- 

c c 


sin 6 - 


X \ cos^ 0 




m 

using the same logic as we used to derive the fractional 
interferometer in section IlYl 

The principles of a generic waveform interferometer are 
essentially the same as those for a regular Fourier inter¬ 
ferometer (Section|n|. The array on the ground measures 
a horizontal slice out of the electric held. The measured 
waveform at a given time is some kind of projection of the 
waveform in time. If the source is inhnitely far away, the 



























15 



FIG. 7: Interferometry can be performed using any set of basis functions, not just Fourier modes. A smaller zenith angle (left) 
stretches out the waveform in space compared to larger zenith angles (middle). Sources that are close effect a “chirp” in the 
waveform, stretching it out more on one side of the array than others (right). The lines styles are the same as in Figure [T] 


projection simply stretches out the waveform depending 
on its altitude. The waveform is stretched out more for a 
source at high altitude (equation [HD) than for a source on 
the horizon (see Section HIl and Figure [7]). By measuring 
the amount the waveform is stretched out on the ground, 
one can reconstruct where on the sky the source is. 

When the source is not at infinity, the generic wave¬ 
form interferometer is analogous to the fractional Fourier 
interferometer (Section ITV| . The source’s position on the 
sky has a measurable parallax from one end of the array 
to the other. Thus, the stretch of the waveform on the 
ground varies from one position to another because of 
the quadratic dependence of the delay time on position 
(equation 15^. On the side of the interferometer closer to 
the source, the waveform is stretched more than on the 
far side. This is illustrated in Figure [7] (right panel). By 
measuring the chirping of the waveform, one gets infor¬ 
mation on its distance. 

As with the usual Fourier interferometer, the time de¬ 
pendence of the electric field signal allows one to dis¬ 
tinguish whether its shape at a given moment is due to 
the radiation’s temporal variation or is some kind of ge¬ 
ometrical effect. We can see from equation |62] that the 
geometry of the array merely shifts the delay time with 
position. Otherwise, the time structure of the signal is 
not changed. 


B. Why the Fourier transform? 

We have greater freedom in interpreting the time struc¬ 
ture of a radio signal than generally appreciated. What 
is so special about the Fourier transform? I believe the 
ubiquity of the Fourier transform arises from four rea¬ 
sons. The first is that the Fourier transform is well stud¬ 
ied, with plenty of “fast” algorithms existing. There are 
widely available free packages that are tuned to com¬ 


pute the FFT quickly. Yet fast algorithms exist for other 
transforms. Fast algorithms now exist for all LCTs, if 
not quite as sophisticated as for Fourier transforms. Ad¬ 
ditionally, the discrete wavelet transform has been known 
for more than 20 years. 

Second, complex Fourier modes are defined as expo¬ 
nentials, which are eigenfunctions of standard derivatives 
and integrals. This simplifies the analysis with standard 
calculus, as opposed to wavelets for example. Yet there 
is no reason, in principle, why we could not use a calculus 
based on different operations. 

Third, the Fourier modes are solutions to the wave 
equations as encountered in physics, like Maxwell’s equa¬ 
tions. Dispersion, and the effective velocity of the signal, 
depends on pure frequency only. Dispersive effects are 
particularly important for radio waves passing through 
interstellar plasma. A mode with several Fourier fre¬ 
quencies disperses. By using the Fourier transform, one 
can account for the effects of dispersion simply. This is 
less of an issue at high frequencies. 

Finally, and I believe most important, the astronomi¬ 
cal sources currently being studied are slowly varying - 
at least compared to one cycle of a radio wave. Radio 
sources have well-defined spectra that can be integrated 
over time. Specifically, the power spectral density of ra¬ 
dio signals from typical astronomical sources does not 
rapidly change in the Fourier domain. As the Fourier 
transform is natural for long-lived sources, there has gen¬ 
erally not been a need for another basis. 

But as the study of transient sources becomes more 
prevalent, there will be greater need for waveforms that 
are well-localized in time (as in Figure [7|). Although the 
Fourier basis can be used for short pulses, the time struc¬ 
ture of the pulse is “hidden” in the phase information. 
The Fourier basis power spectral density of a pulse is 
constant at all frequencies, the same as for generic white 
noise. In fact, the ANITA instrument for measuring radio 
pulses from ultra-high energy neutrinos already employs 
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a kind of “pulse interferometry” for this reason [H^, HH . 


D. The chirpolator as an FrFT 


C. FRBs and the Fractional Fourier Transform 


A possible testing ground for generalized interferome¬ 
ters are FRBs. FRBs are essentially chirps that last for 
a few milliseconds at each frequency; they appear as thin 
curves in “waterfall” plots of the usual time-frequency 
domain M- Assuming an astronomical origin, the dif¬ 
ferent arrival times of the FRBs at different frequencies 
is caused by dispersion in plasma, parameterized by the 
dispersion measure (DM): 



1-1 

Bannister et al. (2011) [8^ noted that the delay time has 
a linear approximation because of Taylor’s theorem: 


, / DM \ ^ V 

^‘”-*‘-’’^“( lOOOpccm-» j(Gife) (Vj- 

(64) 

As I discussed in Section El linear chirps are separa¬ 
ble with a FrFT of the correct order. We must put the 
chirp variable in dimensionless terms to use the FrFT in 
equation 1321 This can be done by using any reference 
frequency vq as the unit for the frequency and l/zzo as 
the unit for time; then the chirping is (1 /VQ){dv/dt). For 
a burst that is smeared out by plasma dispersion, and us¬ 
ing the linear approximation (equation 16411 , the “angle” 
of the FrFT we want is 
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(65) 

For a very small DM, a —>■ 0 because the pulse is basically 
the same at all frequencies. Then we want to examine 
the signal essentially in time domain, where it appears 
as an undispersed pulse. For a very large DM, a —> 7r/2 
because the pulse appears to be a nearly pure tone, which 
very slowly slides in frequencies over time. Then we want 
to examine the signal essentially in frequency domain, 
because the pulse is smeared out to a narrowband tone. 

The difference between a and 7r/2 in eauation l65l seems 
extremely small, but it is in fact measurable. As with 
the fractional interferometer fSection IIVDL the differ¬ 
ence J-a and becomes measurable when the phase dif¬ 
ference between a pure tone and the chirp is ~ 1 radia n. 
The signal must be sampled for a time t ^ 1 / \Jdv/dt Ri 
v-^/^\a--K/2\: 


t>91 fisec 


DM 


1000 pc cm 
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( 66 ) 


Since radio surveys last much longer than a few microsec¬ 
onds, a FrFT is actually quite useful. 


The problem of searching for chirps in signals for a 
sparse array has been considered before by Bannister et 
al. (2011) [1^. Their algorithm for finding chirps in¬ 
volves multiplying the complex voltages at one element 
with the complex conjugate voltage at a second element. 
Since the chirp at one element is just a time-delayed copy 
of the chirp at another element, the product can be shown 
to be a pure tone. A Fourier transform of the product 
finds the frequency of the tone, giving information on the 
chirp. They call this method the “chirpolator”, and note 
that it is the maximum likelihood estimator for chirps as 
found in [s^ . 

Not surprisingly, this operation turns out to be the 
FrFT in disguise. The kernel of a FrFT with a given by 
equation [65] is (c.f., equation (321) 


Ka{t,t') = y/1 — cot a exp[iTTi'Q{{t'^ +t^) cot a+2t't CSC a)], 

(67) 

using the z/q to take care of units. In this case t' has the 
units of time, and we can in fact interpret it as a time 
delay. Thus, the kernel of the FrFT becomes 


Ka{t, —At) = a/ 1 — cot a exp[—i7rz>(t -|- At)^] 

X exp[27r/r'Q(cot a — csca)tAt]. (68) 


Given a chirp with frequency derivative zz, the FrFT mul- 
■ tiplies it by a conjugated, time-delayed copy of itself, and 
then takes a Fourier transform of the result. This is pre¬ 
cisely what the chirpolator does, up to a multiplicative 
constant [s^. A similar algorithm by Bannister et al. 
(2011) is the chimageator, which involves taking the spa¬ 
tial Fourier transform of the voltage product. 

The chirpolator and the FrFT formulations have dif¬ 
ferent strengths. With the chirpolator, it is not necessary 
to know the frequency slide v beforehand to get a reason¬ 
able estimate of the chirp properties. One can apply the 
kernel in equation 1681 without knowing what a is, up to 
a frequency shift. In addition, the antennas do not have 
to be in a rectangular grid to apply the chirpolator. The 
FrFT’s first advantage is that it is conceptually simpler 
- just a rotation in Wigner space (Figure I2). Perhaps 
more importantly, the complexity of a Fast FrFT only 
goes as A^ln for N elements, whereas the chirpolator’s 
complexity goes as N'^. One might consider combining 
these methods’ advantages by using the chirpolator with 
a few elements to get a decent estimate of a, and then 
using the FrFT on all of the complex voltages in the ar¬ 
ray to extract a precise characterization of a chirp with 
all of the data. 
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VII. FULLY GENERALIZED 
INTERFEROMETRY 

A. Interferometry in L2 space 

All of the tools of functional analysis can generalize 
interferometry fully, as they have done for signal pro¬ 
cessing. The electromagnetic field is a complex function 
of position. Formally, we need four real numbers to de¬ 
scribe it, one for each Stokes’ parameter. But for simplic¬ 
ity, just consider a single scalar, the electric held intensity 
at a given point as a function of time, E{t). {E{t) can 
also stand for any other wave, including sound.) Any re¬ 
alistic signal lasts only a hnite time and contains a hnite 
amount of energy: 

/ OO 

\E{t)\'^dt < OO (69) 

-00 

But otherwise, E' may in principle be any square- 
integrable function. 

E is therefore an element of the Hilbert space L 2 , which 
represents all square-integrable functions.® It has an in¬ 
ner product dehned by 

/ OO 

f{t)g{t)dx. (70) 

-OO 

We can dehne various orthonormal bases for L 2 from or¬ 
thogonal functions with {f\g) = 0. The vector of E in 
L 2 in any basis is given by its projection onto each basis 
vector using the inner product. The important part is 
that we can change between bases. 

The standard basis for E is the Dirac position basis, 
with orthonormal vectors representing each 6{t — a). In 
this basis, the components of E are just given by its value 
at each position, as can be seen by plugging 5{t — a) into 
equation 1701 We can say a trivial interferometer merely 
gives E. The Fourier basis is also common, with basis 
vectors of exp(27rwt). The Fourier basis components of 
E in L 2 is just the Fourier transform of the function. A 
normal Fourier interferometer resolves E into this basis, 
inferring where radiation is coming from. These do not 
exhaust the possibilities. All LCTs on E - including 
the FrFT, Laplace transform, and Fresnel transform - 
perform a basis change on L 2 , where the basis vectors 
are the integral kernel. A major theme of this paper has 
been to explore when these other bases are appropriate. 

The LCTs are just the beginning. Other sets of or¬ 
thogonal functions include polynomials, Bessel functions, 
wavelets [s^, and chirplets. In the most general sense, 
an interferometer projects E into any L 2 basis. The only 
question is whether the bases have any useful meaning, 


The space L 2 is perhaps most familiar as the space of all quantum 
wavefunctions. 


and whether the transform is stable both forwards and 
in reverse (it is not in general; lzi,ill). 

An interferometer array measures a cross section of the 
electric field along a line or plane. If the source radiates 
isotropically, then the different antenna locations detect 
the same signal but sample it with different time delays. 
The radiation from a source at infinity arrives as planar 
wavefronts with no curvature, in which case the time de¬ 
lay is exactly proportional to position. If the source is 
closer, the time delay is roughly a second-order function 
of position. Thus, the signal measured by the array at a 
given instant should be a stretched or chirped version of 
the signal in time. An interferometer decomposes E{t) 
into a set of basis functions, and then finds the inner 
product of these basis functions and the signal measured 
in space to determine the direction and distance of the 
source. 

We use Fourier interferometers because the Fourier ba¬ 
sis is the natural one for the time structure of the electric 
field, which in turn is due to the steadiness of typically- 
studied astronomical sources Isection IVl Al) . The electric 
signals from a distant source have a simple representa¬ 
tion in this basis, making it convenient for calculation. 
It is analogous to reference frames in cosmology: there is 
only one reference frame in which the cosmic expansion 
is isotropic at any point, so we usually do calculations 
in this frame. But the laws of physics are the same in 
all reference frames, and nothing stops us from working 
in any frame we wish. Likewise, although Fourier inter¬ 
ferometers are computationally useful, nothing stops us 
from doing interferometry in any other basis of L 2 . 


B. Nonlinear and quantum interferometers 

I have focused on linear transformations of the E sig¬ 
nal, but nonlinear transformations are also possible. For 
a monochromatic input wave, a nonlinear interferometer 
produces a warped output image. This allows us to dis¬ 
tort the Wigner space distribution of a signal in arbitrary 
ways. This might be used to fit signals with quadratic 
dependence of frequency on position (cubic dependence 
of phase on position), including those from sources that 
are very nearby. Likewise, one could also warp the time- 
frequency basis of a signal to fit quadratic chirps. The 
main obstacle nonlinear interferometer, though, is the 
computational cost of calculating a relevant nonlinear 
transformation. 

Quantum optics is an increasingly important field. 
Photon arrival times is already an integral part of in¬ 
tensity interferometry of thermal sources [s^. Quantum 
non-demolition measurements may also prove useful in 
detecting structures below the diffraction limit [o^. Just 
as classical optics can be emulated with a standard dig¬ 
ital computer, it may be possible to emulate arbitrary 
quantum optical systems with a quantum computer. The 
catch is that when the quantum computer produces an 
answer, the wavefunction collapses or decoheres. Since a 
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quantum state cannot be exactly copied, it is lost for¬ 
ever to us and cannot be recovered. Thus, while we 
can run a classically recorded signal through any num¬ 
ber of simulated classical optical systems simultaneously, 
we can simulate at most one quantum optical system. 
The fundamental limitation is that we cannot determine 
which detector a photon passes through. Equivalently, 
the phase of a single photon cannot be measured [9l|. 


C. SETI and generalized interferometers 

To end on an extremely speculative note, an expanded 
consideration of the possible basis functions of a signal 
may have implications for the Search for ExtraTerres- 
trial Intelligence (SETI) [o^. So far, the most common 
assumption has been that a deliberate message would 
appear as an extremely narrowband signal [93|. Ere- 
quency combs are another possibility j94| . as are very 
short pulses [o^. Messerschmitt argued that dis¬ 
crete points within blocks of time-frequency space could 
encode messages with optimal energy efficiency. 

The Fourier basis seems natural to us, so we naturally 
think of tones and pulses. But with SETI we have to ex¬ 
trapolate to entities with potentially very different psy¬ 
chologies and mathematics far beyond our comprehen¬ 
sion. For example, the exponential functions that make 
up the Fourier basis functions are the eigenfunctions of 
the standard differentiation and integration operations. 
But are these the only possible operators that an ETI 
may base their calculus around? If an ETI has studied 
quantum mechanics for a long time, one could imagine 
that the wavefunctions of the quantum harmonic oscilla¬ 
tor may seem a more natural basis set. Or consider that 
wavelets were essentially unknown to the physics com¬ 
munity [s^ when modern radio SETI was first proposed 

93. 

In time-frequency space, basis functions that seem “ob¬ 
vious” to a long-lived ETI may appear very complicated 
to us if they are not using the Fourier transform (Fig¬ 
ure [HI). A more general way to look for ETI signals is to 
search for strong signals that contain both hne-grained 
(or “thin”) features and big, smooth regions in time- 
frequency space. Narrowband signals, combs, pulses, and 
sets of discrete points, are just a few examples of such sig¬ 
nals. Leaving behind smooth ID curves, we could imag¬ 
ine a signal that looks like a 2D pattern in time-frequency 
space, such as a Sierpinski triangle. Known natural phe¬ 
nomena generally produce signals that are smeared out, 
with no hne-grained features, while white noise appears 
uncorrelated in all bases, with no big, smooth regions. 
Even then, if an ETI uses some nonlinear (disconnected, 
even) mixture of time and frequency, the signals may just 
appear as blobs in time-frequency space. Wavelets, after 
all, use a three-dimensional signal representation . In 
addition, looking for mixed hne-grained and smooth fea¬ 
tures in time-frequency space may be a strategy biased 
towards hnding patterns that human perception empha- 



FIG. 8: The most commonly sought signal in radio SETI is 
a narrowband signal, approximating a tone (a). Rapid pulses 
(b) are another sought ETI signal. But if an ETI does not 
use the Fourier basis, the signal may look very complicated 
to us in time-frequency space (c). 


sizes (such as edges in the visual held), which are contin¬ 
gent on human evolution. 

Of course, it becomes impractical to search for all 
possible signals, unless they are distinguished by some 
other way, like being extraordinarily bright signals from 
a Solar-type star. Most of the possible signals look simply 
like white noise. We may be guided by the assumption 
that the ETIs’ resources are limited. The LCTs and other 
transforms like the wavelet transform are advantageous 
in that they can be performed in only linear computa¬ 
tional time. In contrast, a typical linear transformation 
is a matrix multiplication that requires at least quadratic 
computational time. But ultimately, we may just have 
to follow the example of Bell Burnell in discovering pul¬ 
sars [13 , and examine time-frequency space (or any other 
obvious representation that occurs to us) closely for any¬ 
thing out of the ordinary [c.f.. 1^. 

VIII. CONCLUSIONS 

There is nothing inherently special about the Fourier 
basis that requires its use for interferometry. In time do¬ 
main, the signal can be decomposed into any set of basis 
functions. Basis functions that are more wideband but 
localized in time may be more useful when characteriz¬ 
ing very short transients. The position dependence of 
the electric field measured by a planar array is related to 
the signal in time domain for an isotropic source, which 
is stretched out due to projection effects if the source 
is away from the zenith, and is chirped if the source is 
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nearby. The Fourier modes are natural mainly because 
astrophysical sources vary slowly in time. 

Digital interferometers simulate optical systems. The 
standard Fourier transform performed by interferometers 
focuses light coming from infinity. I have described in 
this paper fractional interferometers, which use the Frac¬ 
tional Fourier transform. These focus light from sources 
that are closer than infinity. I have demonstrated that 
there are many sources that are in the regime where a 
fractional interferometer is useful, including atmospheric 
radio transients, artificial satellites, and VLBI observa¬ 
tions of objects within the Solar System. The FrFT with 
orders very close to 1 can be implemented by combining 
the standard Fourier transform and chirp multiplications. 
I have shown that the fractional interferometer is anal¬ 
ogous to the chirpolator [s^. Other Linear Canonical 
Transforms simulate different optical elements. 

Since interferometers these days record signals digi¬ 
tally, it is possible to simulate optics that are not phys¬ 
ically possible. I considered the possibility of digitally 
emulating negative refractive index superlenses. These 
would be useful for constructing images of sources that 
are within one wavelength of the array, which may be 
useful at very low frequencies. These techniques can¬ 
not be applied to sources that are astronomically dis¬ 
tant, because noise overwhelms any remaining signal 


from evanescent waves. This is related to the imprac- 
ticality of using the Laplace transform as a basis for in¬ 
terferometers - the se par ation of exponentials is generally 
an ill-posed problem jT^ . 

My discussion in this paper has been general and qual¬ 
itative. There are many practicalities that need to be 
considered in practice. Most importantly, it is important 
to extend the results to the sparse arrays typical of inter¬ 
ferometers, instead of the dense arrays used in FFTTs. 
What is required is a generalized version of visibilities 
and the van Cittert-Zernike theorem. Secondly, I have 
ignored the sky’s curvature. To accurately implement 
wide-held general interferometers, fractional and other 
extensions of spherical harmonics (99| would be more ap¬ 
propriate, or the sky could be divided into small facets 
[l^ . In addition, I have ignored the polarization of the 
electric held. Calculations of the signal-to-noise ratio of 
observations by general interferometers in the presence 
of noise must be carried out. 
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